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ABSTRACT. 

New BeppoSAX observations of the nearby archetypical starburst galaxies (SBGs) NGC253 and 
M82 are presented. The main observational result is the unambiguous evidence that the hard (2-10 keV) 
component is (mostly) produced in both galaxies by thermal emission from a metal-poor (~ 0.1—0.3 
solar), hot (kT ~ 6- 9 keV) and extended (see companion paper: Cappi et al. 1998) plasma. Possible 
origins of this newly discovered component are briefly discussed. A remarkable similarity with the 
(Milky Way) Galactic Ridge's X-ray emission suggests, nevertheless, a common physical mechanism. 

1. Main Interests 

The main interest of studying SBGs in the local universe is that it is believed that the 
starburst phase was very frequent during the early (z>l) evolution of galaxies (Madau et 
al. 1996). Thus understanding the properties of local SBGs means, presumably, under- 
standing the properties of most galaxies of the early universe. Studies of SBGs can also 
help to shed light on key astrophysical issues like: i) the AGN/starburst connections, 
since both phenomena require concentration of material to the nuclear regions and are 
possibly linked to merging phenomena and/or presence of molecular bars (Maiolino et 
al. 1998), ii) the chemical enrichment of the intergalactic medium (IGM), because of the 
frequent presence in SBGs of chemically-enriched SB-driven outflows (Heckman 1997) 
and iii) the structure, evolution and even formation of galaxies through the effect of the 
"feedback" provided by massive stars (Ikeuchi & Norman 1991). 

2. Current Observational Status 

Because of the short space available, we limit the following discussion to the 2-10 keV 
energy band (see Fabbiano 1989 for a review on the optical and soft X-ray properties of 
SBGs). In this energy range, the best-quality data available to date have been essentially 
limited to the two brightest SBGs NGC 253 and M82 (but see recent studies on lower- 
luminosity galaxies by Delia Ceca et al. 1997). Data from the Einstein IPC, EXOSAT 
ME, Ginga, BBXRT and ASCA (Ptak et al. 1997 and ref. therein; Tsuru et al. 1997, 
Moran & Lehnert 1997) have clearly demonstrated the existence of a hard component 
well described by either a thermal model (kT ~ 6-9 keV) or a power-law model (T ~ 
1.8-2.0). The absence in the data of a significant Fe-K line emission has, however, always 
been puzzling and has induced several authors to propose explanations alternative to 



the thermal emission, i.e. the presence of a LLAGN (Ptak et al. 1997, Tsuru et al. 1997) 
or inverse-Compton (IC) scattered emission from the interaction of IR photons with 
SNe-generated relativistic e~ (Moran & Lehnert 1997). 

In the following, we present new BeppoSAX results for these two SBGs which demon- 
strate that such hard component is most probably emitted by a thermal plasma. In this 
paper, spectral results showing the first evidences of Fe-K line emission (at ~ 6.7 keV) 
and high-energy steepening expected in the case of thermal emission are presented for 
both NGC 253 and M82 (see Persic et al. 1998 for more details on the results for NGC 
253). In a companion paper (Cappi et al. 1998), we combine these spectral results to 
the (unprecedented) spatial information of BeppoSAX between 2-10 keV. 

3. BeppoSAX spectra of NGC 253 and M82 

BeppoSAX observed NGC 253 on Nov.29-Dec.2, 1996 for - 50, 113 and 50 Ksec and 
M82 on Dec.06-07, 1997 for - 29, 85 and 30 Ksec, for the LECS, MECS and PDS 
operating between 0.1-4 keV, 1.3-10 keV and 13-60 keV, respectively. A Galactic column 
density of 1.28 x 10 20 cm" 2 for NGC 253 and 4.27 x 10 20 cm" 2 for M82 were added 
in all the following spectral fits (Dickey & Lockman 1990). The spectral analysis was 
performed using version 10.00 of the XSPEC program. 
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Several continuum models (e.g. a power-law, a bremsstrahlung, a Raymond-Smith, a 
multi-temperature Raymond- Smith, a Mewe-Kaastra-Liedahl model) and combination 
of them with and without absorption, with solar and variable abundances have been 
used to model the spectra (Fig. 1). We have found that both spectra are best described 
in terms of a two-component thermal model (Fig. 2). At energies lower than ~ 3 keV, 
emission lines (e.g. FeL, SiXV and SXV) characteristic of a soft thermal component 
arc detected, in agreement with ASCA results. At higher energies, both sources require 
the continuum emission to steepen (thus excluding a power-law description) and, most 
significantly, an Fe-K emission line at E ~ 6.7 keV is detected in both objects (see Table 
1 and Fig. 1). Interestingly, both hard components need to be absorbed in order not to 
over-produce the continuum at E <1 keV (as also found with ASCA, Ptak et al. 1997). 

Best-fit results obtained by fitting both source spectra with a two-temperature Mewe- 
Kaastra-Liedahl model with variable abundances ("vmekal" in XSPEC) are shown in 
Table 2 and Fig. 2. Following Persic et al. (1998), in order to reduce the number of 
free parameters, the abundances of He, C, and N were fixed at the solar value and the 
heavier elements divided into two groups: Fe and Ni (most likely associated to SNc I 
products), and the a-elements O, Ne, Mg, Si, S, Ar and Ca (most likely associated 
to SNe II products); elements in the same group were constrained to have common 
abundances in solar units. It should be noted that although some extra contribution 
from a non-thermal (hard) component cannot be formaly excluded, we estimate that it 
should not contribute more than <~ 20% of the total 2-10 keV flux. 



The detection of the Fe K line allows us to reliably determine the Fe abundances 
of the line-emitting matter: we find a value of ~ 0.25 and 0.06 solar for NGC 253 and 
M82, respectively, consistent with the sub-solar values (based on upper- limits) predicted 
from previous hard X-ray observations (see §2). Although more insecure because of the 
current uncertainties of theoretical models in computing the Fe-L emissivity, similar Fe 
abundances are found for the soft component. On the other side, a-elements seem to 
have larger abundances (~ 1.7 solar for NGC 253 and ~ 2.4 solar for M82). 



Table 2: Two-components (vmekal) thermal models 
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4. On the origin of the hard, thermal, component 

On one side, our results for the soft X-ray component arc in good agreement with 
previous X-ray results on SBGs (see §2), i.e. that it originates from a SB-driven galactic 
superwind (see Persic et al. 1998 for more details). On the other side, the origin of the 
hard thermal component is more puzzling. As demonstrated in Cappi et al. (1998), it 
is also extended possibly flowing out of the galaxy plane (as for the soft X-rays) in the 
case of NGC 253. Among possible candidates, one could consider a collection of point 
sources (e.g. X-ray binaries, SNRs, clusters of protostars) or truly diffuse emission due 
to IC scattering of IR-optical photons from relativistic e~ or, finally, a very hot (~ 
6-10 x 10 7 K) ISM phase. Alternatives to the hot ISM phase, however, could hardly be 
reconciled with the observed spectrum because i) the average spectra of X-ray binaries 
would predict a power-law continuum and Fe K line at 6.4 keV, ii) the average spectra of 
SNRs would predict too softer (kT< 4 keV) spectra, hi) clusters of protostars have too 
low luminosities (L 2 -iokcV ~ 10 32 erg s _1 ) and iv) IC emission would predict a power-law 
spectrum. The most natural explanation seems therefore a hot thermal diffuse emission, 
i.e. a "super super-wind" most likely powered by the starburst. 

It is interesting to note the striking analogy of the present results with the ASCA 
observations of the Galactic Ridge (GR) (compare our Fig. 2 with Fig. 3 in Kaneda et 
al. 1997): also in this case, hot diffuse (and absorbed) emission is detected at a similar 
temperature and, if scaled for a ~ 1/100 times lower Galactic star formation rate, 
with similar energetics. The only difference being in the somewhat larger (~ 0.8 solar) 
metal abundances deduced for the GR, these remarkable similarities strongly suggest a 
common physical origin. 

5. Conclusions 

The main conclusion from the above results is that the BeppoSAX observations of 
the two archetypical SBGs NGC 253 and M82 have revealed for the first time strong 
evidence of a hot, really hot, thermal diffuse plasma responsible for their hard (2-10 keV) 



0.5 1 2 

channel energy fkeV) 



12 5 10 

channel energy (keV) 



Fig. 1. Spectra of NGC 253 and M82 fitted with a single power-law model (r ~ 1.7) to illustrate 
the spectral complexities (emission lines and continuum curvature) 
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Fig. 2. Best- fit spectra of NGC 253 (left) and M 82 (right) with a 2T Mekal model (Table 2). 



emission, and that this discovery may have profound implications on AGN/starburst 
connections, on the chemical enrichment of the IGM and on the formation and evolution 
of galaxies. 
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